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Abstract 

We prove that Mathieu's N — 2 supersymmetric Korteweg-de Vries 
equations with a = 1 or 4 admit Hirota's n-supersoliton solutions, whose 
nonlinear interaction does not produce any phase shifts. For initial profiles 
that can not be distinguished from a one-soliton solution at times t <^ 0, 
we reveal the possibility of a spontaneous decay and, within a finite time, 
transformation into a solitonic solution with a different wave number. This 
paradoxal effect is realized by the completely integrable N = 2 super-KdV 
systems, whenever the initial soliton is loaded with other solitons that are 
virtual and become manifest through the r-function as the time grows. 
Key words and phrases: Hirota's solitons, N — 2 supersymmetric 
KdV, Krasil'shchik-Kersten system, phase shift, spontaneous decay 

Introduction. Let u{x, t; 9i, 62) be a scalar complex bosonic N — 2 superfield. 
We consider Mathieu's N ~ 2 supersymmetric Korteweg-de Vries equations [TJ 

m 

d d 

ut = -u^xx + 3(uPiI?2u), + ^{ViV2u')^ + 3au'w,, V,^— + e,- — 

(1) 

with a = 4 or a = 1. In this paper, we construct n-supersoliton solutions of ([T]), 
M = A-I?i2?2logTfc,,...,fc„(a;,i;6'i,6l2), ^(a) = const e M, (2) 

that possess the following properties. First, the linear truncations of Hirota's 
exponential series, 

n 

r = 1 + aj exp rjj , rji = kiX+Lu(ki)-t + X{ki)-0i92, /c^ G M, G R, (3) 

1=1 

yield exact solutions of ([T]). Functions ^ contain no pairwise interaction 
{Aij(ki^ kj) = 0, see ^ below) and no multiple-order interaction terms (re- 
spectively, Ai^...i^{h^, fcj,) = ^ic'hi^ic^hp) = 0, see [3]). Second, 
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as n different solitons constitute an n-soliton solution and obey the nonlinear 
superposition formulas while overlapping, they acquire no phase shifts having 
become distant from each other. Thus the spatially extensive solitary waves 
demonstrate a dimensionless behaviour of particles with elastic nonlinear col- 
lisions. We recall that a similar example of the field representation for point 
particles was known, e.g., for soliton solutions of the integro-differential Ben- 
jamin-Ono equation [3]. 

Third, we show that a fast soliton that precedes a set of other waves always 
decays spontaneously and, sooner or later, the observed solution converts into 
the slowest soliton from that set. We say that the initial soliton, which was 
visible at t ^ 0, was loaded with n — 1 virtual solitons ahead of it. We demon- 
strate how the fastest soliton gains the virtual ones and becomes virtual itself, 
while the slowest virtual soliton manifests itself behind the others and becomes 
visible. Of course, these properties of the new solutions for = 2 super-KdV 
equations ([T]) are radically different from superposition laws and phase shifts 
that are exhibited by the n-solitons for the standard KdV equation, see ([6|) be- 
low, and by the super-soliton solutions ([7]) derived from them. This effect, which 
is new to the best of our knowledge, may be useful in verifying the relevance of 
nonlinear N = 2 supersymmetric KdV equations ([T]) in field theory and, more 
generally, establishing the physical sense of the Grassmann variables Oi, 92- 

The paper is organized as follows. First, we summarize well-known proper- 
ties of n-soliton solutions for the KdV equation. Here we also recall how these 
solutions are inherited by Hamiltonian N = 2 super-KdV systems ([T|) from the 
purely bosonic KdV. Then we formulate our main result, describing the r-func- 
tions and dispersion relations in the new n-supersolitons for ([T]). If a = 1, then 
the new super-fields yield exact n-soliton solutions of the Krasil'shchik-Kersten 
equation [5]. We establish the paradoxal properties of the ri-solitons for ([T]), 
including their spontaneous decay followed by transition into virtual states. 

§1. Known n-solitons for the N — 2 super-KdV. Some classes of exact 
travelling wave solutions u{x + ct+i6i92, 0i+i92) for equations ^ with a — —2, 
±1, ±4 were obtained in [Sj using the symmetry invariance technique. First, a 
Lie subgroup of point symmetries for ([T]) with a G M was chosen such that the 
respective group invariants specified the arguments of the solutions u as above. 
It follows easily that the Taylor expansion u = po + {9i + 16*2) • pi in the second 
argument, with both po and pi depending on x + ct + i9i92, converts ((T]) to the 
triangular system of ODE upon po and pi . The second assumption of [3] is that 
solutions of the first equation in this system, which is pg — apQ + i{a + 2)poPo + 
cpo + const ~ 0, have no movable critical points. This yields the admissible 
values a — —2, ±1, ±4. Particular solutions po{x + ct + i9i92), pi = 0, which are 
found in Eq. (3.15)] for ([T]) with a e {1,4}, are close to the one-solitons ^ 
with c = u>/k in However, the pairwise or collective interaction of the 

travelling waves was not discussed in [S] . 

At the same time, it is clear that the N = 2 super-KdV equations ([T]), and 
not only with the values a € {—2, 1,4} that ensure the complete integrability, 
possess multi-soliton solutions. Indeed, expand the bosonic N ~ 2 super-field u 
in 9i and 6*2, 

u = uq + 9iui + 92U2 + 9i92Ui2, (4) 
and substitute ([4]) in ([1]). Then, in particular, the even components uq and U12 
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satisfy the system (TJ [21 [S] 

uo-t + uq-xxx - Sauluo-x + (a + 2)(uoUi2)x - (a - 1)(uiM2)x = 0, (5a) 

+ 6ui2Ui2-x — (a + '2)uoUO:xxx — 3a(uo;2:UO:xa; + UqUi2-^x + 2lioMO;xUl2) 
- 3u2M2;a:K - (« + 2)uiUi-xx + 6a{uoUiU2-x ~ UoUi.xU2 + U0;xUlU2) = 0. 

(5b) 

Set the odd components ui = U2 = 0. First, if a = —2, then system ([S]) 
spHts to the modified KdV equation (|5ap upon uq together with the KdV-type 
equation (j5bp with a feed-back. In what foUows, we do not consider the well- 
studied case a = —2. If, under the same assumption ui = U2 = 0, we let a = 1, 
then system ^ is the Krasil'shchik-Kersten coupled KdV-mKdV equation [S], 
see also [3 El E] . In the next section we describe a new class of its multi-soliton 
solutions in Hirota's form. 

Conversely, set the components mq, wi, and U2 to zero. Then for any a e R 
the entire N = 2 supersymmetric equation ([1]) reduces to the purely bosonic, 
scalar KdV equation 

Ul2-t + U12-XXX + Gui2Ui2-x ^ 0- (6) 

Therefore the super-function 

U = 0102 ■ Ui2{x,t), (7) 

where ui2 satisfies (jSj), is a solution of ([T]). Reciprocally, each function ([7]) yields 
the solution ui2 — 'D2'Di{u) of ([6|). Note that for the class ^ one transfers the 
formulas of Backlund autotransformation for KdV to the N = 2 super-KdV 
equations (|T|). 

Substitution ([7]) further implies that equations ([T]) inherit the classical Hirota 
n-solitons dD) from the solutions 

ui2= A{a) -(f/dx^ log T (8) 

of the KdV equation fB]), see [3], where the r- function is 

n 

T = ^ CXp(^^ m ■ {r]i + logai) + ^ fliflrn ■ log AijTL^ , (9) 

fj.e=0,l 1=1 l<l<m<7i 

l<f.<n 

and the dispersion and interaction (now A = 0) are given through 

A{a) = 2, Lu = -k^, Aera = ih-k^f/{h + kmf, aeeR. (10) 

Lemma. After each pairwise interaction, which is encoded by Agmiki, km) 
in PH)) . the solitons for KdV equation ^ acquire the phase shifts, whose abso- 
lute values depend on the wave numbers ki, km and equal 

Aifim ^ 2log(^\ke - km\/\ke + km\y, (11) 

the signs of the shifts themselves are opposite for the ^-th and m-th solitons. 
For instance, if ki > k„i > 0, then the fast soliton with the wave number ki is 
pushed forward, and the slow soliton given by fc,„ retards back. As t —^ —oo, the 
two solitons move freely, still staying (respectively, as t +oo, having already 
become) distant from each other. 
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The knowledge of this effect dates back to the beginning of the KdV-boom in 
late 60's, see [3] and references therein. We are obliged to re-derive formula pT|) 
in order to compare it with the application of the same reasoning to the new 
solutions (frf]) . see below, of = 2 super-equations ([I]). The new super-solitons 
demonstrate a principally different behaviour both from the analytic and phys- 
ical points of view. 

Proof. Formula (jlip is particularly transparent on the level of the r- functions ([9|) ; 
without loss of generality, we assume n = 2 and thus let ^ = 1, m = 2 for a 
two-soliton solution of © given by ki > k2 > 0. Let xi{t) and X2{t) be the 
coordinates of the peaks of the first and second solitons, respectively, at the 
times t when they are sufficiently distant from each other. For convenience, we 
set all ai :— 1, and hence in view of ([9]) the r- function has the form 

'rkik2{x,t) = 1 + exp?7i + cxpry2 + ^12(^1, ^2) • exp(r/i -I- r/2). 

First, let i ^ be fixed. Consider the fast first soliton, which has not yet 
overpassed the slow second one. In the vicinity of its peak located at xi{t), 
we have 771 w and, since the peak of the second soliton is still far ahead at 
2^2 (i) 3> xi{t) for i ^ 0, we deduce 772(2;) ^ for a; « xi{t). Hence0 

Tfcife2(x,t) = H-expryi +expr72 • (l + yli2 -expTyi) =t I + exp?]! =rfei(?7i), 

(12) 

because exp 772(2;) JJ, uniformly as t — s- ~oo on any finite interval that con- 
tains xi(t). The asymptote is just the fast soliton in its pure form ^ without 
any phase shift yet. 

Next, consider a vicinity {x ~ 2:2 (i)} of the peak of the slow soliton, which 
still goes far ahead the fast one. Now 772 ~ and, since X2{t) » xi{t) for t <^0, 
we have rji (x) ^ at any x ^ X2 (t) . Consequently, 

Tfeifc2(x,t) = exp?7i • (^1 + exp(r]2 + logAi2) + exp(-r/i) • (l + exp7y2' 

Acting by the operator 2d^/dx^ olog and noting that the phase rji is only linear 
in X (it does not contain higher powers of x, we conclude that 

Ufcifc,(x,t) = 2— ^(fcia;+wii)+log(^l+exp(772+2 1ogf^)-|-exp(-77i)-(l+exp7/2 

since exp(— 771) Jj. on a finite interval around X2{t) as t —00. The negative 
phase shift is precisely as t ^ +00, it will vanish, and hence the overall 

phase shift of the slow soliton will be positive. 

Notice that, initially, the slow soliton is displaced forward by jAi/j^j^jl with 
respect to an identical soliton ([5]) with the same wave number /c2 that would 
arrive to the origin X2{t) = at t = 0. Thus, we recall, the system of coordinates 
(x, t) used in Hirota's method is such that the center of mass for the two soliton 
solution achieves 2; = at some t < 0. 



^ Here and everywhere below the standard symbol =t denotes the uniform convergence of a 
functional sequence on a finite closed interval. The symbol JJ, denotes the uniform convergence 
from above (in our case, to the zero function). 
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Second, suppose t +00; by now, the fast soUton has left the slow one far 
behind. Therefore, in a vicinity {x w X2(t)} of the slow peak with 'q2{x) ~ 0, 
we have rjilx) ^ at a fixed time t, whence 

Tfcifc2 {x, i) = 1 + exp r/2 + exp r/i • (l + A12 cxp 772) =^ 1 + cxp r/2 = r^^ (m), 

because exp7yi(a;) H near X2{t) as t ^ +00. 

Concentrated around xi{t) such that 771 = 0, the fast soliton shifts forward, 
borrowing the energy and momentum from the slowly moving obstruction. In- 
deed, by r]2{xi{t)) we have 

Tk^k2{x,t) = exp 772 ■ (1 + exp(?7i +logAi2) + exp(-?72) • (l + exp77i)^, 

which implies 

d^ r / 
Uk,k2{x,t) = 2— |^(fc2X+a;2i)+log(^l+cxp(77i+21og|^^)+exp(-7;2)-(l+exp7/i^ 

^Ufe,(77i+21og|^) (13) 

on {x « xi{t)} as t ^ +00, because exp(— r/2) Ji there. Comparing the in- 
and out-going states of the two solitons, we recognize the accumulated phase 
shifts ±Aipkik2 for the first and second solitons, respectively. □ 

Remark 1. The dispersion law in (|10p is determined by the Hirota bilinear 
equation for the KdV equation see j^, 

D:,{Dt + Dl)T-T = 0, 

where and Dt are Hirota's derivatives. We recall that a bilinear represen- 
tation for the Manin-Radul N — 2 supersymmetric KdV equation ([T]) with 
a = —2, see [10], was obtained in [TT] for its expression in Inami-Kanno's 
form [12]. The bilinear representation for the a = —2 super-KdV (|T]) is such 
that the t- functions for the n-soliton solutions, assigned to the wave numbers 
ki, . . ., fc„, contain the Grassmann constants Ci, ■ • Cn- 

§2. New (virtual) rt-solitons for the N — 2 super-KdV. We discard the 
idea of using any 'Grassmann constants,' see Remark [T] above. Consequently, 
the most general admissible form of the Hirota exponentials is 

r] ^ kx + Lut + Xei02, (14) 

where £ K are the wave numbers, uj are the frequencies, and A S C de- 
scribes the dependence on the Grassmann variables 6*1, 02 for n-supersolitons 

Uki...k„{x,t;di,92). 

Proposition 1. i. The travelling wave solutions with Hirota's exponen- 
tials for the N = 2 supersymmetric KdV equation ^ with a = 1 
or o = 4 are 

u = A • PiX'2log(l + aexp??), (15) 

where the dispersion and the relation between the wave number k and the 
nonzero factor A are 

Lu = -fc3, A = ±ik. fc e R, a e R; 
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ii. The Hirota two-soliton solutions of ^ with a G {1,4} and A 7^ in p^ . 



u = A - T>i'D2 log(l + ai exp?7i + 02 exp 772 + 74i2aia2 e^pivi + ^2)), 

exist if and only if relations ([1^1 for Hirota's exponentials hold and 
there is no coupling between the exponentials for one-soliton solutions (fT^ : 
Ai2(fci,A:2) =0. 

The exact n-supersoliton solution of the N ^ 2 super-KdV ^ with a G 
{1,4} and A, 7^ is 



1 + 2^i=i o-i exp -qi dx 



(17) 



Here n G N is arbitrary and we put the reduced phase be 

f\ = kx + uit; (18) 

then, substituting (fTH)) for A(a), w(A:), and A(fc), we expand the supersoli- 
ton ^ with the r-function 

En 
Ui expryj • (1 + \i9i92) 
i—i 

in 01^2, whose square is zero. 

iv. The a — —2 super-equation ([1]) does not possess travelling wave so- 
lutions (fT51) other than PH)) with A = 0, and hence it does not ad- 
mit n-supersolitons (jl7p for any n > 1. This also holds for any other 
a G K \ {1, 4}, which still do admit the n-solitons given by ([7HTU|l but not 

by dUHlel). 

Proof. Once ([TMT7| are known, statements HHml are verified by direct calcu- 
lation. The exhaustive conclusion of statement [iv] is obtained using analytic 
environment [13j . □ 

Remark 2. In [T3], a particular class of exact solutions was constructed for 
the Manin-Radul N — 2 super-KP equation and for its reduction to the N = 2 
super-KdV equation ^ with a = — 2. The components of these solutions (those 
are called solitinos in loc. cit.) are at most quadratic in Hirota's exponentials 
exprji for any number n > 2 of interacting solitinos, see [T4, Eqs (39-40)]. 

We have extended this result onto the entire triple a G {—2, 1, 4} of integrable 
Hamiltonian super-equations ([1]) . Namely, for N — 2 super-KdV equations with 
a = 1 or 4 we revealed the n-soliton solutions p7p . whose r- functions remain 
linear in exprji for any n>l. We stress that the unexpected ri-solitons p7p of 
the N = 2 super-KdV must also solve the N = 2 supersymmetric KP equations 
that reduce to ^ with a G {1,4}. Definitely, the admissible super-systems are 
not the Manin-Radul N = 2 super-KP [TD] that implies a = — 2. 

Finally, we note that the constraint ui = U2 = is fulfilled for the compo- 
nents of super-fields (fT7|) . Therefore Proposition [T] yields exact solutions of the 
Krasil'shchik-Kersten equation, see (O, whenever a = 1. Likewise, it provides 
the n-solitons for the bosonic limit of ([T]) if a = 4. 
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Under the nonlinear superposition of n solitons, exact solutions (|17p of super- 
equation which are given by linear truncations ^ of Hirota's t- functions, 
demonstrate the properties that are essentially different from the standard pic- 
ture (fTTj) for KdV (HI), regarding both the asymptotes at the time infinity and 
the phase shifts for each of the n solitons. We claim that 



the asymptotic behaviour (t ±00) of an n-soliton solution (fT7|) 
for ^ is completely determined by the soliton that moves behind 
the others, and that no phase shifts at all are accumulated as the 
time t runs from —00 to +00. 

The proof goes in parallel with (fTMOj) . see p. HI and the reasonings are alike for 
all n > 2, hence we set n = 2 and suppose fcmax = ki > k2 = fcmin > 0. Without 
loss of generality, we assume ai = 1 for all i. Recall that Tkj^k2ix,t] 61,62) — 
l-|-exp77i-|-exp 772, where the full phases rji are and the dispersion is specified 
by (|16p . This yields the solutions pT)) of where A{a) is given in ([TC]) for 
either a = 1 or a = 4. 

Notice that the essential non-constant part of a soliton (flSl) is, as usual, 
concentrated near the zero of the reduced phase (fT8|) . However, the first, purely 
imaginary summand in the expansion (|17p does not have a peak located at 
?7 w 0. At the same time, the real coefficient of 6162 in the one-soliton solution 
does have the peak at 77 « 0. (Let us remark that this coefficient is nothing 
else but one half (for a = 1) or a quarter (respectively, a = 4) of the one-soliton 
solution ([8]) for KdV, see ([T0|) .) We shall use this in what follows to refer on 
points X G M. 

So, let t —00. In a vicinity of the peak xi{t) of the fast soliton with 
the maximal wave number ki — fcmax, the reduced phase IjlSp is ryi « 0, and 
772 ^ 0. For still distant solitons, exp 772 JJ, as i — > ^00 on any finite interval 
{x xi{t)}. Therefore, 

Tk^kA^,t]6i,e2) = l-fexpryi -|-exp772 =4 l-l-exp??! =rfe^(77i). (19) 

This is the fast soliton with its unshifted full phase rii(x,t). 

On the other hand, near the peak X2 (t) of the slow soliton we have 772 ~ 
and 771 3> 0, whence 

Tkik2ix,t;6i,62) = expiji ■ (l -I- exp(-?7i) • (1 + exp?72)). 

Applying the operator A{a) ■ Pil?2 ° log to the r-function, we obtain the asymp- 
tote of the two-soliton solution on a finite interval near X2{t): 



Uk,k2 {x, t- 61,62)^ A V1V2 \^{kix+LJit) +log(l+Ai6li6l2) +log (^l+exp(-77i) ■ (l-t-exp 772) 

=t —A\i = const. 

In other words, the solution near X2(t) does not vary under a small finite motion 
left or right along x at t Q. For this reason, we say that the slow soliton has 
not yet manifested its presence, remaining virtual. We also say that the visible 
fast soliton, which starts behing the slower soliton(s) in the beginning of the 
elastic scattering process, is loaded with the virtual soliton(s). Note that the 
asymptotic behaviour of Uk-i---k„ix,t]6i,62) at early times t <C is analogous 
for any n > 2. 
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Now consider t — > +00. The slow soliton with the wave number /c2 = fcmin 
becomes visible. Near its peak at X2{t), we have 772 and 771 <C as soon as 
the fast soliton has gone sufficiently far ahead. Consequently, 

Tk^k2{x,t;6i,92) = l + exp?7i +exp?72 =^ l + exp?72 ^Tk^im), 

because exp?7i JJ, near X2{t) as i ^ +00. 

And the fast soliton becomes virtual. Indeed, near its peak at xi{t) such 
that fji « 0, we have 7)2 ^ and exp(— 772) j]. as t ^ +00, whence 

Tk^k2{x,t;0i,92) = exp 772 ■ (l + exp(-?72) • (1 + exp77i)) 

yields 

Ukik2ix,t;0i,92) = AVi'D2 {k2X+LU2t)+log{l+X20i02)+log(l+cicp{~r]2)-{l+e-Kprii 

=1 -AX2 = const. (20) 

This value is the limit lim^^+00 Uk2 (2;, t) determined by the slow soliton that 
now goes behind the fast. The constant does not alter under a small finite 
motion left or right from the point 771 = 0, where the peak of the fast soliton 
would be located at any t :s> 0- 

We have described the "vanishing" of the fast soliton. Now the visible slow 
soliton is loaded with the virtual fast soliton. 

It remains to note from (fT^H^n]) that none of the solitons, neither the visible 
fast that becomes virtual, nor the virtual slower soliton that becomes visible, 
acquires any phase shift. The proof is complete. □ 

Let us summarize our main result. 

Theorem 2. i. As t — > — cxd, the time asymptote of the n-soliton (I17p for 
the super -KdV ([1]) is given by the fastest soliton with the maximal wave 
number fcmax, which is loaded with the virtual slower solitons that corre- 
spond to k2, . . ., fcmin < fcmax- These virtual solitons are invisible att <^ 0. 

ii. As t +00, the asymptote of the same solution Uk-^...k„{x, t; 9i, 02) for ([T]) 
is the slowest soliton with the minimal wave number fcmin, which is loaded 
with faster solitons determined by fcmax, ■ • fcn-i > fcmin- These n — 1 vir- 
tual solitons remain {for fc2, . . ., fcn-i) or become (respectively, for fcmax) 
invisible at t ^ 0. 

Hi. None of the solitons, neither the visible, which goes behind the others, nor 
the virtual ones in front of it, acquires any phase shift under the evolution 
in time t £ (—00, +cxd). 

Corollary. The iterated spatial and temporal limits of the pure imaginary first 
summand in (|17p are not permutable: as x,t +cx), one haS 

lim Hm exp77,^ ^ A (fcmax) 7^ A (fcmin) lim hm exp77,^ 

t^+oox~,+oo I + ^^^^aiexprji a:^+oo^ f^+oo 1 + ^^^^^ exp77i 



^Note that the limit in x in the right-hand side is taken before the peak x„{t) of the slowest 
n-th soliton passes through a point x, which means that the exponentially vanishing tail of 
all solitons has not arrived yet. 
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This also shows that the solution Uki---k„ix,t:6i, 62) is not a Schwarz function 
in X. Clearly, as x — > —00 at any t, the limit of the first summand in (jl7p equals 
zero (we assume ki > everywhere in the text). The real coefficient of 9i92 
in (fT7|) tends to A{a) d/dj:(const) = as a; — > ±00 at any t. 

The initial limit as a; — > +00 at any finite time t is due to the visible fastest 
soliton with the wave number fcmax, which is loaded with virtual slower soli- 
ton(s). The spatial limit as a; ^ +00 for the out-going state at large times 
< ^ is specified by the slowest soliton, which is loaded with the virtual faster 
waves and which progressively occupies the entire a:-axis. 

Remark 3. The a = 4 super-KdV equation ^ possesses twice as many Hamil- 
tonian functionals and hence twice as many symmetries as the other two cases 
a G {— 2,f}. Therefore there are twice as many commuting flows as in the 
hierarchies with a = —2 or a = f . In particular, between the spatial transla- 
tion d/dx and the translation d/dt, which is determined by ([T]) in the a — 4 
super-KdV hierarchy, there is the N — 2 supersymmetric 'Burgers' fiow [f 5 , 

ii — 'Di'D2Ux + uux- (21) 

Proposition [T] and Theorem [2] are literally reproduced for (|2ip . The n-superso- 
litons are p?)) with the following dispersion relations in Hirota's r- function ([3]), 

A = 2, uj = -ik^, X = ik, a G M. (22) 

Conclusion. The new n-soliton solutions (fT7)) satisfy N = 2 super-KdV equa- 
tion ([1} with a G {1,4} or, to be even more precise, its bosonic limit, which 
is the Krasil'shchik-Kersten system if a = 1. Recall that the t —00 asymp- 
tote of our solutions (fT7|l is given by the fastest soliton with the maximal wave 
number fcmax, which is loaded with n — 1 virtual slower solitons ahead of it. 
We have proved that, looking at the initial profile and recognizing it as a one- 
soliton solution up to any precision of measurements, one can predict neither 
the out-going state nor even its asymptote —iAkmin at t as a: — > -l-oo. 

For an observer, the fast soliton (fci = fcmax) is subject to a spontaneous 
decay. The profile transforms into a collection of nonlinear-overlapping waves 
{ki > k2 > ■ ■ ■ > kn) that, finally, constitute what one sees as another soliton 
with a different wave number fc„ = fcmin- For a given number of solitons n > 2 
and for a priori given wave numbers fcmin and fcmax, this can be realized along 
any trajectory, which is a point in the configuration space 

R!;:-^ XM"-1 3 (fc2,fc2 -fc3,---,fcn-2 -fcn-l) X (a2,...,a„), 

since over C one can always rescale ai = 1. 

The absence of the phase shifts, supplemented with the presence of hierar- 
chies of conservation laws for Hamiltonian N = 2 supersymmetric KdV equa- 
tions ([I]), demonstrates that the spatially extensive solitary waves behave in 
elastic collisions as dimensionless particles (as massive material points) . In this 
paper, we established the possibility of their spontaneous decay, and we formu- 
lated the laws for transformation of observed particles into the virtual ones and 
vice versa. 
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